Path-integral Monte-Carlo simulations for electronic dynamics on molecular chains: 

II. Transport across impurities 
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Electron transfer (ET) across molecular chains including an impurity is studied based on a recently 
improved real-time path integral Monte Carlo (PIMC) approach [J. Chem. Phys. 121, 12696 
(2004)]. The reduced electronic dynamics is studied for various bridge lengths and defect site 
energies. By determining intersite hopping rates from PIMC simulations up to moderate times, the 
relaxation process in the extreme long time limit is captured within a sequential transfer model. The 
total transfer rate is extracted and shown to be enhanced for certain defect site energies. Further, 
it is revealed that the entire bridge compound approaches a steady state on a much shorter time 
scale than that related to the total transfer which allows for a simplified description of ET along 
donor-bridge-acceptor systems in the long time range. 



I. INTRODUCTION 



Charge transfer in native or designed molecular struc- 
tures has attracted a substantial amount of research re- 
cently This is mostly due to the perspective 
to exploit corresponding processes in a controlled way in 
integrated nano-electrical circuits 0, Q . However, even 
though tremendous progress has been achieved, still var- 
ious rather fundamental questions need to be elucidated. 
Certainly, one of these is the role of the environmen- 
tal degrees of freedom, originating from solvent and/or 
residual vibronic molecular degrees of freedom, which are 
often neglected in approaches based on transport the- 
ories for clean mesoscopic devices. Nevertheless, they 
are known to be crucial for charge transfer in condensed 
phase structures 0, H, ■ This is assumed to be partic- 
ularly true in cases where impurities along the transport 
channel are to be surmounted by an intimate interplay 
between charge-phonon interactions such as thermal ac- 
tivation, nuclear tunneling, and coherent electronic tun- 
neling (superexchange) . 

The development of effective approaches to include 
complex environments in the quantum dynamics of elec- 
tronic systems has been one of the major goals in chem- 
istry and physics as well . The widely used and pow- 
erful Redfield formulation has revealed important aspects 
of charge transfer in various systems 0, 0, 0, 
However, as an approximate method it is limited to weak 
friction with high frequency vibronic modes and suffi- 
ciently high temperatures. In a previous paper [T5l |. re- 
ferred to as Ref. [I] henceforth, we used a numerically ex- 
act Monte Carlo path integral approach (PIMC) to cap- 
ture electronic transfer (ET) dynamics in donor-bridge- 
acceptor (DBA) systems. In essence, it exploits that 
the path integral formulation allows for an exact elim- 
ination of the bath as well as the quasi-classical elec- 
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tronic degrees of freedom which reduces the dimension- 
ality of the Hilbert space considerably (T^ . The method 
has the merit of being applicable in the whole parame- 
ter space, that is also in those ranges where approximate 
methods fail, e.g. for low frequency environments and low 
temperatures. It further allows to treat bridges with a 
larger number of electronic sites and up to sufficiently 
long times so that in Ref. [I] we could investigate the 
length dependence of the total transfer rate from donor 
to acceptor, which has also been the subject of recent 
experiments J 7,, ,18] ■ 

In Ref. [I] we paid particular attention to charge dif- 
fusion across a DBA complex with energetically degen- 
erated bridge states. Here, we continue this study by 
analysing the impact of an impurity located in the cen- 
ter of the bridge. One may expect that the existence of a 
defect is the typical situation in native degenerated struc- 
tures. Moreover, meanwhile advanced chemical synthesis 
allow to "dope" molecular chains in a controlled way and 
study transport properties by linking them to external 
leads |23|. The archetypical case is thereby that of a 
defect with an energy gap to its surrounding electronic 
sites as this scenario includes effectively also defects in 
the electronic intersite coupling. Our main focus here lies 
on a microscopic description of the transport dynamics 
and the total transfer rate, which for isolated DBA sys- 
tems corresponds to the relaxation rate towards thermal 
equilibrium, and for integrated DBA structures subject 
to a dc- voltage is related to the steady state flux . In 
general, the long time limit where a global thermal equi- 
librium state is established, cannot be reached within 
typical times accessible in PIMC simulations, however, 
our improved method covers time windows from which 
all dynamical scales including the relaxation rate can be 
extracted. Thus, we use the exact PIMC data to deter- 
mine and verify a sequential transfer model and, based on 
this procedure, obtain a complete understanding of the 
dynamical features of incoherent ET in DBA structures. 

The paper is organized as follows. We briefly out- 
line the model for describing dissipative ET dynamics 
in Sec. and in Sec. IIIII the corresponding real-time 
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FIG. 1: Molecular chain with 6 = 5 bridge states {d — 7) 
and an impurity in its center. Gray bars denote the different 
energetic layouts for the impurity investigated in this paper. 



PIMC simulation approach; for further details we re- 
fer to our previous paper Ref. [I]. Then, the numeri- 
cal results are presented in Sec. IIVI where we start in 
Sec. IIV iCI bv extracting intersite hopping rates from the 
PIMC data based on a sequential hopping model. The 
long time dynamics according to these findings is anal- 
ysed in Sec. IIV Bl and the dynamics of the total bridge 
compound in Sec. IIV CI At the end a short summary and 
some conclusions are given. 



II. DYNAMICS OF THE DISSIPATIVE d-LEVEL 
SYSTEM 

We investigate the dynamics of a single electron mov- 
ing on an impurity flawed molecular chain. The gen- 
eral approach to describe single charge transfer across 
bridges has been outlined in Ref. [I]. Here, we thus col- 
lect only the key ingredients. The molecular chain con- 
sists of d = 25* + 1 discrete sites (see Fig. 0, where, 
separated by equal distances a, b = d — 2 adjacent bridge 
states {l-Bl), . . . , \Bb)} are clamped between a donor \D) 
and an acceptor state \A). The latter two are energeti- 
cally degenerated and define the 0-energy level. The bare 
bridge states are also energetically degenerated but at a 
higher energy hesi = = 2.5?iA, while the impurity, 
located in the center of the bridge |/) = |-B(d+i)/2), ex- 
hibits an energy gap e/ — ^B^J^^l)/2 ~ with respect to 
the other bridge states (we only regard chains symmet- 
ric with respect to the impurity site, i.e. with an odd 
number of bridge sites). The focus of our study lies then 
on the transport across the bridge/impurity compound 
when this gap is energetically raised or lowered. 

We restrict our model to nearest-neighbor coupling, 
that is, electronic motion is facilitated through tunneling 
which only occurs between adjacent states with a uniform 
real tunneling amplitude l^s,s' ■ The electronic coordinate 
can then be expressed as 



where —S < s{t) < S. The position operator thus is 
equivalent to the spin S operator 



aS^ |s} = as\s) 



(2) 



with \s) denoting the (orthonormal) localized electronic 
states. In terms of electron transfer, \—S) and \S) rep- 
resent the donor and acceptor \D) and \A), respectively, 
while the other states are referred to as the bridge states. 

The free d-level system (dLS) Hamiltonian can then be 
written as Ref. [I]. 



(3) 



where describes the energetic distribution of the elec- 
tronic sites according to 

E,|Z?) = E,|A) = 0, 

E,|S,) = eslB,) for 1 < z < 6 , {d+l)/2 , 
E,|/) = {eB+€i)\I) , (4) 

while Sx governs the tunneling, 

^x\D) = A|Bi) , 
S,|Si) = ^{\D) + \B2)) , 

Sx|S(d+l)/2Tl) = ^(|5(d+l)/2T2) + > 
S,|/) = A(|B(d+i)/2_i) 

+ l^(£;+i)/2+i)) , 
Sx\B,) = A(|B,_i) + |B,+i)) 

for i e {2,. ..((i+ l)/2 - 1, 
(d+l)/2 + l,...,d-l} , 
Sx\Bb} = AQA) + \Bb-i)) , 
Sx\A) = A\Bb) , 

(5) 

The interaction with solvent and vibronic degrees of free- 
dom is introduced within the framework of a systern-plus- 
reservoir model, leading to the total Hamiltonian 



H — i?dLS + Hj + He 



E 



p2 

2nia 



(6) 



Here, the residual degrees of freedom are archetypically 
modeled as a harmonic bath coupling bilinearly to the 
position of the electron (Hi) As discussed in detail 
m Refs. [ISISlEl, 

this provides a reasonably accurate 
description of reality for the great majority of ET sys- 
tems. 

For the dissipative electronic dynamics we focus on the 

time- dependent site populations, 



q{t) = a ■ s{t) , 



(1) 



P,,,,,(t) = Tr{e* 



Ht/hi 



/)(*■/ 1 



(7) 
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which are normaHzed, ^s/,si (i) = I7 and where 

Wi(0) specifies the initial state of the total compound. In 
most theoretical and experimental works an initial sep- 
aration of the electron and the environment is assumed 
corresponding to an initial density matrix 



(8) 



with the electron held fixed in state \ si) and the bath nor- 
malization Zb = Trje^'^-^^} assuring the full system's 
density matrix to be normalized for all times. For a trans- 
fer process across the entire chain one typically prepares 
the electron initially in the donor state, i.e. \si) = \D) 
Above, n is the dipole moment of the electron, 
and £ denotes the dynamical polarization of the bath 
[IC/I, which is equilibrated with respect to the initial po- 
sition of the electron. By comparing with Eq. lIBll. we see 
that /Lt£ — aY^^CaXa- As pointed out in Ref. [2^, this 
"standard preparation" often used in ET experiments is 
especially suitable for a theoretical description of thermal 
transfer rates. 

Switching to the path integral representation of 
Eq. O, the bath degrees of freedom can be traced out 
exactly, yielding the reduced dynamics 



P. 



At) 



T^s Ss{t),sf exp i tScIlsIs] - 



(9) 

Here the path integration runs over closed paths s{t) 
starting at s(0) — Si and propagating along the real-time 
contour t e ^ t ^ 0, which connect the forward and 
backward paths s{t') and s'{t'), respectively. Further- 
more, S'(jLs['S] denotes the total action of the free dLS. 
The influence of the traced-out bath is completely en- 
coded in the Fcynman-Vcrnon inRuence functional $[s] 



$[5,5']= dt' dt"[s{t')- s'{t')][L{t' -t")s{t") 
Jo Jo 

- L*{t' ~t")s'{t")] 

+ ,|^V[.2(t')-,s'2(0], (10) 

which is written in terms of the complex- valued bath au- 
tocorrelation function reading for real time t 

1 r-,,^(,)CoshNy/2 .t)] 



sinh(?i/3w/2) 



with /3 = l/ksT and 



il= - doj—^ = 2Aci 

'0 



TT 



(12) 



Both quantities are completely determined by the spec- 
tral density 



(13) 



which for a condensed-phase environment typically be- 
comes a continuous function of In the classical regime, 
the description of the bath reduces to the celebrated Mar- 
cus theory '3, ^ , where the dissipative influence onto the 
electronic dynamics are fully governed by the classical 
reorganization energy HAd. 



III. SIMULATION METHOD 

We perform path-integral Monte Carlo (PIMC) simu- 
lations, utilizing a discretized version of the path-integral 
expression (0 where the electronic system's path s{t) is 
replaced by {sj}i<j<2q, Sj = s({j — l)t/q). The "dynami- 
cal sign problem" [23 notorious for real-time calculations 
of quantum systems is significantly relieved by employ- 
ing an optimized filter based on the blocking approach 
which exploits certain symmetries of the influence func- 
tional 16] . Besides some improvements described below, 
we used the same computer code as described in detail 
in Ref. [I]. 

The computational bottleneck of this approach is the 
creation of the MC trajectory since the evaluation of each 
trial move requires the calculation of a series of matrix 
multiplications. However, here significant saving can be 
achieved. Starting with the standard MC weight. 



WMc[{£.j},Vq+i) = exp 



E 

k>j=-. 



(i^W(6,...,C,)---i^^'He,)) 



2si,r(,+ i 



(14) 



where the exponential is the real part of the influence 
functional while the second factor represents a product 
of complex matrices, which include the complex valued 
contributions from the influence functional and the prop- 
agation of the undamped electronic system (for details, 
see Ref. [I] App. A), in a first step, these matrices are 
stripped off all dissipative terms leading to 



exp 
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E 

k>j=2 



^k^k-jCj 



(15) 



i^(o,6)---^(e„o) 



Now the matrices K{£^j, can be computed from the 
free propagator alone and thus are, unlike the X*^^^ ma- 
trices in Eq. p4(l . Jocai in time in the sense that they 
connect only two neighboring time slices. Accordingly, 
the former can be calculated and stored beforehand for 
all combinations of £,j and £,j+i, while the latter have to 
be recalculated every time a trial move is evaluated. 

The respective speed up can be multiplied by comput- 
ing and storing sequences of K matrices prior to perform- 
ing the MC moves, e.g. with 



K(")(6,...,en+l) 



n^(^„o+i), (16) 
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which diminishes the number of matrix multipUcations 
necessary to obtain wmc by a factor of n. Of course this 
approach is Umited by its memory consumption which 
increases exponentially with d. However, with the on- 
board memory facilities of nowadays computer's being in 
the GByte range, it still offers tremendous performance 
improvements in terms of simulation duration. 

Finally, in a last step the X^"' are replaced by real- 
valued matrices which are obtained from the former by 
replacing all matrix entries with their modulus. Similarly 
to ignoring the dissipative terms present in the original 
K^^") matrices, this somewhat impairs the statistics of the 
PIMC simulation. However, in both cases the thus nec- 
essarily larger number of accumulation processes is more 
than compensated by the speed up of the evaluation of 
the trial moves, provided the system is not too far from 
the coherent-incoherent transition. Fortunately, far away 
from this transition the application of special techniques 
to soothe the sign problem usually is not necessary any- 
way. Note also that the deviation of the improved MC 
weight H15|) from Eq. H14|l decreases with increasing n. 

To illustrate the effectiveness of this approach we per- 
formed calculations both with the improved as well as 
with the original MC weight. For the parameters used 
in this paper, the former achieved a speed up of typi- 
cally a factor of 12; the data for e//A = 5 presented in 
Fig. 13 e.g., were obtained on IBM-P615 within a CPU 
time of 6.2 days with n = A (cf. Eq. ((TB|l and q — 100 
discretization steps on the real-time axis, where 6.7 x 10^ 
accumulation processes were performed; the statistical 
error in the average sign was 1.2%. 



?i/3A = 0.1, which gives still rise to strong quantum ef- 
fects in the bath {uich(3 = 0.5) though. 

Following these lines, we expect an initial nonequi- 
librium state to decay - after an initial transient on a 
time scale rtrans ~ multi-exponentially towards thermal 
equilibrium. According to our previous simulations [I], 
the equilibrium occupation probabilities are taken to be 
Boltzmann distributed, i.e.. 
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(18) 



Here, Pg',P^,P^, andP^, i = 1, . . . {d+l)/2-l,{d + 
l)/2 -I- 1, . . . , 6 denote the populations of donor, acceptor, 
impurity, and the remaining bridge sites, respectively, for 
t oo. For intermediate to long times this decay is 
governed by a mono-exponential relaxation with a time 
constant that is the inverse of the total transfer rate Ft- 
The crucial question is then, to what extent signatures 
of the impurity can be detected in Ft, which often is 
the only observable accessible in experiments. In case 
that the molecular chain is connected to external leads 
being integrated into an electrical circuit, an applied dc- 
voltage gives rise to a steady state with a flux that is 
directly related to Ft 0- 



IV. TRANSFER DYNAMICS 

For vanishing dissipation, the charge dynamics of an 
electron initially localized at the donor site displays co- 
herent oscillations corresponding to no net-transfer to the 
acceptor site. These oscillations survive for weak cou- 
pling to an environment and sufhciently low tempera- 
tures, but with decreasing amplitude as the system re- 
laxes towards its equilibrium state. In the sequel, we 
are interested in completely incoherent transport, where, 
however, quantum mechanical properties of the phonon 
bath play a substantial role and give rise e.g. to nuclear 
tunneling. According to our previous study, we consider 
a spectral density of the form 

J(w) = 27racje-"/"'% (17) 

where the damping strength is chosen as a = 1/4, well 
below the coherent-incoherent transition a ~ 1/2, and 
the cut-off frequency Wc/A = 5. As seen in Ref. [I], for 
these parameters coherent oscillations in the electronic 
dynamics show up for h(3A > 0.3. Therefore we fix the 
inverse temperature sufhciently above this transition at 



A. Population dynamics and intersite hopping 
rates 

For molecular chains with energetically degenerated 
bridge states, but variable length treated in Ref. [I], a se- 
quential hopping model turned out to capture the popu- 
lation dynamics for times beyond the transient time scale 
quite well. Further, direct transfer populating bridge 
sites only virtually, known as superexchange, turned out 
to play no substantial role for the bath parameters spec- 
ified above even for larger energy gaps between the in- 
termediate and donor/acceptor sites. Consequently, we 
assume first a purely sequential model to be applicable 
also here which is self-consistently verified by analysing 
the QMC data. The influence of superexchange is in- 
cluded in a second step. 

Then, the population dynamics for P = 
(Pd , Pbi , • ■ ■ , , -Pa) is governed by 

P^AP, (19) 

with a rate matrix independent of time for t > rtians 
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A 







-^BD — 
Tb 





Tb 

-2rB 





Tb 



\ 



Tb 




-^B — Tb/ 




Tib 
-2r/B 



Tb/ 

-^Bi — 







-2Ti 

Tb 




Tb 

—^B — Tba r^B 

Tba -r^B / 
M 



Note that Eq. H19|l can be derived from an exact non-local 
retarded Master equation ^2&]- Due to the sequential 
hopping, the intersite rates Tdb, Tbd = {Pb^ /Pd)^db, 
and Tb should be independent of the gap energy hej of 
the impurity, as our simulations confirm (see below) . The 
latter one determines only Tbi and F/b which are related 
by the detailed balance condition 

pco 

Ts/ = -T^^iB = e^'^P'TiB . (21) 

The intersite rates were then obtained from our QMC 
data by the following procedure: For increasing values 
of the transient time scale rtrans starting at Ttrans = 
, the populations stemming from solving Eq. H19(l with 
integration constants specified by the QMC results for 
-Pa (Ttrans), . . . , Pd (n rans) wcrc fitted to the QMC data 
for PA(t > Ttrs^ns), ■ ■ ■ , Poit > Ttrans)- This fixed the 
intersite rates as functions of the transient time scale. 
Each true rate constant then was extracted as the plateau 
value of its corresponding rate function with rtrans as the 
actual transient time. The rates were fitted in sequel 
rather than at once. First, Fbib and Fb were obtained 
from a fit to a completely degenerated bridge, i.e. e/ = 0, 
with Tif3A = 0.1 and eb = 2.5 A, yielding 

Tdb/A = 0.298, 
Fbd/A = 0.383, 
Fb/A = 0.370, (22) 

which basically coincide with the golden rule values, see 
Ref. [I], Eq. (22). For impurity bridges, e/ ^ 0, the 
bridge- impurity rate Tbi remains thus as the only pa- 
rameter in the sequential population dynamics (|19|l . A 
subsequent fit yields Fb/ as a function of the impurity 
offset e/. The existence of time independent intersite 
rates (plateau time) and the accuracy of this procedure 
proofs in turn the applicability of the hopping model. 

We start by presenting as an example in Figs. [3 El 
PIMC data for the population dynamics of a bridge with 




FIG. 2: Electronic populations for a — 1/4, u>c/A = 5, 
ft/3A = 0.1, es/A = 2.5 and ej/A = 0,2.5,5,7.5, 10, and 12.5 
(circles, squares, diamonds, triangles up, triangles down, and 
triangles right, respectively). 
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FIG. 3: Electronic populations for a — 1/4, ujc/A = 5, 
hl3A = 0.1, es/A = 2.5 and ei/A = 0,-2.5,-5,-7.5,-10, 
and —12.5 (circles, squares, diamonds, triangles up, triangles 
down, and triangles right, respectively). 



d — 7. Our simulations included bridges with 6 = 3,5,7 
(d = 5,7,9) and gap energies between —12.5 < e//A < 
12.5 in steps of 2.5. To obtain rate constants with suf- 
ficient accuracy, the simulations had to be performed 
with exceptional small statistical error and up to long 
times At = 20. For example, for a typical coupling 
A = 300 cm~^ this corresponds to simulation times up 
to 60 ps. In Fig. ^ we display corresponding results for 
Tbi as a function of the energy gap e/. Qualitatively, 
the numerical rates follow the golden rule formula, but 
deviate, as expected, from the classical Marcus expres- 
sion [see Ref. [I], Eq. (21)]. Quantitatively, however, they 
are smaller than the golden-rule results by about 10% or 
less, a fact that can partly be attributed to adiabatic 
effects like recrossings in the Landau- Zener region since 
our bath is not strictly in the nonadiabatic limit where 
the golden rule expression is valid. Rates for different 
bridge lengths differ by less than 2% outside the max- 
imum and by about 5% at the maximum where those 
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FIG. 4: Bridge-impurity rate Vbi as a function of the impu- 
rity energy gap e/. Rates extracted from PIMC simulations 
(markers) for various chain length with d = 5, 7, 9 are shown 
together with the results according to the classical Marcus 
(dashed) and the quantum mechanical golden rule (solid) ex- 
pressions. 



for shorter bridges exceed those of longer ones. In the 
classical limit, the pronounced maximum of VBii^i) cor- 
responds to an activationless transfer at e/ — —Ad, while 
for UfiiOc > 1 collective quantum effects in the bath pro- 
mote nuclear tunneling between adiabatic surfaces, shift- 
ing its location towards smaller gap energies (cf. fig. EI . 
Note that Vbi > F^ around the maximum. Deviations 
between rates for different b are the largest close to this 
maximum where adiabatic effects are expected to be the 
strongest as the bath is directly located around the tran- 
sition point. 

The non-monotonic behavior of FB/(e/) is already well- 
known from Marcus theory since the activation barrier 
between bridge and impurity state is a quadratic func- 
tion of e/ -I- A, see Ref. [I], Eq. (21). Therefore, Marcus 
rates show perfect symmetry with respect to e/ = —A ~ 
— 2.5A. The fact that acceptor populations in Fig.|2and 
121 however, exhibit a symmetric behavior with respect to 
ei — 0, simply reflects the symmetry properties of the 
total rate Fs/(e/) + F/B(e/) = FB/(e/) + TBii-ej). For 
long times, of course this symmetry in the -Pa(0 broken 
due to the different equilibrium populations. 

Based on these results and Eq. (|19|) . we show in Fig. [51 
populations Poit) and Pi{t) for & = 5 at fixed times but 
for varying e/. This gives deeper insight in how the final 
equilibrium populations (|18|l are established during the 
time evolution. Namely, the impurity population tends 
to establish an exponential dependence oc exp(— ?i/3e/) on 
a much shorter time scale than the donor one displays a 
behavior cx exp[?i/3(eB -I- e/)]. We will see in detail below 
that this is a characteristic facet of the DBA dynamics. 
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FIG. 5: Donor population Pd (left) and impurity population 
Pi (right) for d = 7 at fixed times vs. e/. Left: At = 5, 
(solid), 10 (long-dashed), 20 (short-dashed), 100 (dotted); 
right: At = 5, (solid), 10 (long-dashed), 20 (short-dashed), 
50 (dotted). 
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FIG. 6: Total transfer rate Ft according to the intersite 
hopping rates extracted from the PIMC simulations and the 
rate dynamics Eq. (1191 for molecular chains with d = 5, 7, 9 
sites vs. the impurity gap ej. 



B. Total transfer rate 

With TjB at hand, the transfer dynamics according 
to Eq. H19|) is known for arbitrary times. In particu- 
lar, the least non- vanishing eigenvalue of the rate matrix 
A defines the total transfer rate Tt being typically, e.g. 
for longer bridges, much smaller than the intersite rates. 
This way, we get access to information in the extreme 
long time range by extracting intersite rates from exact 
PIMC simulations. 

Results are depicted in Fig. El for varying gap energy 
£/ and bridge lengths 6 = 3,5,7. As discussed, the over- 
all values of Tt are smaller by about a factor of 5 to 10 
compared to the intersite hopping rates meaning that the 
time scales for local dynamics are well separated from 
global relaxation. Strikingly, there is a maximum at 
e//A « —1.5 which is a signature of the activationless 
regime found in the Tgi rate. Hence, a defect with an 
appropriate negative gap enhances the global transport 
dynamics compared to the situation where it is absent. 
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FIG. 7: Total transfer rates without (Ft, dashed) and in- 
cluding superexchange across the impurity site (Fts, solid) 
vs. ei. 

Of course, for very large negative and positive gaps the 
chain is basically cut into pieces and the sequential trans- 
fer to the acceptor tends to vanish. 

In our previous studies, non-sequential electronic 
transport processes (superexchange) played only a mi- 
nor role. For the range of gap energies considered in the 
PIMC simulations, the same holds true here. Since the 
rates gained from the PIMC data are close to the golden 
rule rates, particularly for higher/lower impurity sites, we 
analyse the situation for extreme gap energies by addi- 
tionally including direct transfer without populating the 
impurity, i.e. from \B(^a+i)/2-i) to |-B(d+i)/2-Hi), into the 
rate matrix A via a superexchange rate Tse- Knowing 
that superexchange is relevant for high lying impurities 
only, we estimate TgE from a quantum mechanical golden 
rule expression as well (forth order perturbation theory 
in A), see Eqs. (23), (25) in Rcf. [I]. The corresponding 
^Tsi^i) is shown together with F^ in Fig. and verifies 
that for the total transfer, sequential hopping dominates 
by far up to impurities with e//A < 20. Superexchange 
dominates for larger gaps though, where the transport to 
the acceptor is almost frozen. In this latter regime elec- 
tronic tunneling through the impurity exceeds the effect 
of nuclear tunneling in the sequential hopping mecha- 
nism. 



C. Bridge relaxation 

Even though the molecular bridge including the de- 
fect plays a key role in the charge transfer from donor to 
acceptor, typically in experiments individual bridge site 
populations are not accessible. Thus, the question arises 
whether the bridge compound considered as a sub-unit of 
the entire molecular DBA structure exhibits a character- 
istic dynamics on its own. For that purpose, we look at 
the total bridge population Pb = where the sum 
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FIG. 8: Time dependent populations of bridge Psit), donor 
Poit), and acceptor PA{t) for a bridge with 6 = 5 sites and 
an impurity gap e//A = 5. Dots represent PIMC data, solid 
lines the full population dynamics according to Eq. I19L and 
grey lines the asymptotic behavior according to Eqs. 1231 and 
1241 1. see text. The dotted lines indicate the equilibrium pop- 
ulations. 




FIG. 9: Time dependent populations of bridge Psit), donor 
Poit), and acceptor PAit) for a bridge with 6 = 5 sites and an 
impurity gap e// A = 5 where the impurity is located one site 
to the left of the center of the bridge. Dotted lines indicate 
the respective equilibrium values. 



runs over all bridge sites including the impurity. Clearly, 
due to normalization Pb = 1 — Pd — Pa ■ 

Upon closer inspection one finds that due to the sym- 
metry of the rate matrix in Eq. (|19|l . every second term 
in Psit) drops out exactly so that its asymptotic dynam- 
ics is governed by the second least eigenvalue of A, de- 
noted by Ttb, instead of Tt- For times beyond a certain 
threshold tb ^ I/Ft, the bridge saturation can thus be 
captured by 



PB{t) =P^ + [Pb{tb) - P^] e-i^--(*-^«) 



(23) 



where Txb ^ Tj^; accordingly, even if the intersite rates 
in A are not known, Ttb can still be found by study- 



ing the saturation limit of the time dependent expression 
TTBit) = PB{t)/\PB{t) - Pb\- This theoretical predic- 
tion from the hopping model is completely verified by the 
simulations as shown in Fig. |S1 where Pb (t) for fixed gap 
6/ /A — 5 is depicted as a function of time. Apparently, 
it tends to saturate on a considerably smaller time scale 
than the total chain relaxes towards thermal equilibrium. 
An explicit value for the transient time tb is gained from 
the condition that the actual PB{t) and its asymptotics 
differ by less than 1%. This leads to Atb = 6.4 
which is considerably below the simulation time. As seen 
in Fig. IIUI the bridge relaxation rate takes its maximum 
not at negative, but rather positive gap energies, where 
the overall transfer is slowed down. In this case, the de- 
fect facilitates as a scattering center first relaxation to 
its left side and then to its right one by blocking the flux 
to the acceptor site. Note again that the asymptotic dy- 
namics (I23II of the total compound does not mean that 
the individual bridge site populations have approached 
their equilibrium values, too. 

Now, for very long times t ^ I/Fts one has Poit) — 
—PA{t) so that then the transport across the DBA com- 
plex can be described within an effective two state model 
with a transfer rate being identical to the total transfer 
rate Tt extracted above, i.e.. 



Pd]' 



ro/A) (24) 



Pnit) - Pd+[Pd{td/a) 
PA{t) = P^ ^[Pd{td/a)^Pd\^- 

Here, besides Tt only information about the donor en- 
ters. Further, Ar^i/^ = 17.4 (extracted similar as de- 
scribed above for r;,) obeys tu/a ^ 1/Fts and lies also 
well within the simulation time window, thus verifying 
that the PIMC simulations cover the entire relevant time 
range. Everything beyond this range is in a sense "sim- 
ple" and completely determined by the information re- 
trieved from the numerical data. The comparison be- 
tween Figs. 1^ and [TUl reveals that indeed Ttb ^ Tt- 

While this strong time separation is no longer true for 
a defect located off center, the total bridge population 
comes still significantly closer to its equilibrium value 
on a time scale much below 1/Tt than Po and Pa, see 
Fig. El The remaining deviation falls off with Tt, but 
practically Eq. still gives an accurate description of 
the D-A transfer at long times. Hence, the underlying 
scenario for the DBA transfer is this: On an interme- 
diate time scale 1/Ttb a local steady state for the to- 
tal bridge population is established, which then acts for 
larger times merely as an electronic reservoir mediating 
between D and A. In turn, the comparison between the 
dynamics of PA{t)+P]:){t) and the individual populations 
Pd (t) , Pa {t) would allow experimentally to gain informa- 
tion about the symmetry of the bridge/defect structure. 
In this context it is worthwhile to mention that a sym- 
metric situation treated here numerically is not just an 
idealized toy model but can actually be realized in ex- 
periments, see Ref. poj . 

The above findings on the time scale separation allow 
for a simple recipe to successively extract rate constants 
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ei/A 

FIG. 10: Bridge relaxation rate Ttb vs. e/ for bridges with 
6 = 3, 5, 7 (from top to bottom), see text for details. The solid 
lines are based on the PIMC data, while for |e//A| > 12.5 the 
dashed lines are extrapolations based on golden rule hopping 
rates to better see the fall-off. 

from the donor and acceptor dynamics only. Namely, in 
a first step Tt is obtained from the asymptotic behav- 
ior of Poit) a-iid PA{t)- Then, for a symmetric system 
Poit) + PA{t) allows to obtain the second least rate con- 
stant Ttb- Further, by subtracting the dynamics corre- 
sponding to these rates from Poit) and PA{t) the third 
least eigenvalue shows up etc. Practically, this procedure 
works for the rate constants governing the dynamics from 
moderate to long times for which typically a sufficient 
separation of time scales is guaranteed. 

V. SUMMARY AND CONCLUSIONS 

Based on a recently improved real-time QMC scheme 
to treat numerically exact electronic transfer dynamics in 
condensed phase, we analysed the impact of an impurity 
on transport across molecular chains. Experimentally, 
this is not only the typical situation in native structures, 
but advanced molecular synthesis even allows to assemble 
corresponding chains in a controlled way. We restricted 
our study to the case of an impurity site with an energy 
gap to its neighboring sites, which can be seen as the 
archetypical scenario as it provides e.g. also an effective 



description of a defect in the electronic intersite coupling. 

Our main focus was the total transfer rate correspond- 
ing to the time scale for global relaxation towards thermal 
equilibrium. Since for longer bridges this time scale is in 
general considerably larger than the largest times acces- 
sible in real-time PIMC simulations, we used the latter 
ones to extract, beyond a transient time, intersite hop- 
ping rates between the electronic sites. This information 
was then used in a sequential transfer model, verified by 
the PIMC data, to capture the extreme long time dy- 
namics. 

Roughly, two dynamical ranges can be distinguished in 
the charge transfer: In a first step, the total bridge pop- 
ulation approaches (symmetric bridge/defect) or comes 
close to (asymmetric bridge/defect) its value for thermal 
equilibrium facilitated by the presence of an impurity, 
while in a second step overall relaxation sets in associated 
with a characteristic unique time scale, the total trans- 
fer rate. Even though the total bridge population tends 
to saturate in the first range, intersite relaxation be- 
tween individual sites still takes place. Nevertheless, the 
bridge compound including the impurity can basically 
be treated as a thermal reservoir mediating the transfer 
from donor to acceptor. The extent to which a time scale 
separation between bridge and total relaxation can be ob- 
served is directly related to the bridge/defect symmetry. 
The total transfer rate exhibits a pronounced maximum 
at negative values of the impurity energy gap, thus re- 
flecting the enhancement of the local bridge-impurity rate 
in or near the activationless regime. In principle, a con- 
trolled doping with a defect site may therefore accelerate 
the transfer across a degenerated chain. As in our pre- 
vious study, superexchange across the impurity plays a 
dominant role for extremely large gaps where the trans- 
port comes almost to rest. Our analysis may give qual- 
itative and quantitative insight into experiments similar 
to the one reported in Ref. 20] where a molecular chain 
contacted with leads was designed to include a defect at 
its center. 



Acknowledgments 

We acknowledge financial support from the DEC 
through Grant No. AN336-1 and the Landesstiftung 
Baden- Wiirttemberg gGmbH. J. A. is a Heisenberg fellow 
of the DEC. 



[1] J. Jortner and M.A. Ratner (eds.). Molecular Electronics, 

(Blackwell Sci., Oxford, 1997). 
[2] J. Jortner and M. Bixon, eds.. Adv. Chem. Phys. 106, 

107 (1999). 

[3] C. Dekker and M.A. Ratner, Phys. World 14, 29 (2001). 
[4] A. Nitzan, Ann. Rev. Phys. Chem. 52, 681 (2001). 
[5] P. Hanggi and M.A. Ratner (eds.), Chem. Phys. 281, 
(2002). 



[6] C. Kergueris, J. -P. Bourgoin, S. Palacin, D. Esteve, C. 
Urbina, M. Magoga, and C. Joachim, Phys. Rev. B 59, 
12505 (1999). 

[7] J. Reichert, R. Ochs, H.B. Weber, M. Mayor, and H.v. 

Lohneysen, Appl. Phys. Lett. 82, 4137 (2003). 
[8] R.A. Marcus, J. Chem. Phys. 24, 966 (1956). 
[9] R.A. Marcus and N. Sutin, Biochim. Biophys. Acta 811, 

265 (1985). 



10 



[10] U. Weiss, Quantum Dissipativc Systems, Scries in Mod- 
ern Condensed Matter Physics, Vol. 2 (World Scientific, 
Singapore, 1998). 

[11] A.K. Felts, W.T. Pollard, and R.A. Friesner, J. Phys. 
Chem. 99, 2929 (1995). 

[12] A. Okada, V. Chernyak, and S. Mukamel, J. Phys. Chem. 
A 102, 1241 (1997). 

[13] W.B. Davis, M.R. Wasielewski, M. Ratncr, V. Mujica, 
and A. Nitzan, J. Phys. Chem. 101, 6158(1997). 

[14] D. Segal, A.Nitzan, W.B. Davis, M.R. Wasielewski, and 
M.A. Ratner, J. Phys. Chem. 104, 3817 (2000). 

[15] L. Miihlbacher, J. Ankerhold, and C. Eschar, J. Chem. 
Phys. 121, 12696 (2004). 

[16] C.H. Mak and R. Egger, Adv. Chem. Phys. 93,39 (1996). 

[17] M. Magoga and C. Joachim, Phys. Rev. B 56, 4722 
(1997). 

[18] W.B. Davis, W.A. Svec, M.A. Ratner, and M.R. 

Wasielewski, Nature 396, 60 (1998). 
[19] B. Giese, J. Amaudrut, A.-K. Kohler, M. Spormann, and 



S. Wessley, Nature 412, 318 (2001). 

[20] M. Mayor, C. von Hanisch, H.B. Weber, J. reicliert, D. 
Bcckmann, Angew. Chemie 41, 1183 (2002). 

[21] D. Chandler, in Liquids, Freezing and the Glass Transi- 
tion, ed. by D. Levesque, J. P. Hansen, and J. Zinn-Justin 
(Elsevier Science, North Holland, 1991), Les Houches 51, 
Part 1. 

[22] X. Song and A.A. Stuchebrukhov, J. Chem. Phys. 99, 
969 (1993). 

[23] L. Miihlbacher and R. Egger, J. Chem. Phys. 118, 179 
(2003). 

[24] R.P. Feynman and F.L. Vernon, Ann. Phys. (N.Y.) 24, 
118 (1963). 

[25] Quantum Monte Carlo Methods in Condensed Matter 
Physics, M. Suzuki (ed.) (World Scientific, Singapore, 
1993), and references therein. 

[26] R. Egger, C.H. Mak, and U. Weiss, Phys. Rev. E 50, 
(655) (R) (1994). 



